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ICARUS AND THE VARIATION OF PARAMETERS 
By SAMUEL HERRICK 


1. The minor planet 1566, Icarus, also known 
as 1949 MA and as Baade’s object, possesses 
special interest because of (a) the solar parallax 
and related quantities, (b) the mass of Mercury, 
(c) the relativity effect. Of all minor planets 
known it passes closest to the sun and Mercury. 
Of all minor planets still “‘in captivity”’ it passes 
closest to the earth. And because of its high 
eccentricity it presents a severe test both for the 
relativity motion of its perihelion and for certain 
aspects of perturbation determination, e.g., for 
the method of variation of parameters that is to 
be illustrated herein. 

(a) For the solar parallax and the mass of the 
earth-moon system, whether determined trigono- 
metrically or by the more likely gravitational 
method, Icarus in 1968 will pass the earth almost 
at its minimum distance, about four million miles. 
Thus, roughly, its geocentric parallax will be four 
times that of Eros in 1931 and the gravitational 
effects will be sixteen times as great. Balanced 
against these favorable factors we find Icarus 
generally too faint to be observed except for short 
intervals, and then, unless it is very close, only 
by the largest telescopes. It is often badly placed 
for observation. In 1951 it was too low in the sky 
after the end of astronomical twilight for obser- 
vations to be obtained in either the northern or 
southern hemisphere. In 1952 it was too low in 
the northern hemisphere, but happily five obser- 
vations were obtained at Pretoria, so that the 
orbit of Icarus is now well established. In the 
moderately close approaches of 1953, 1958, and 
1959 Icarus will again be better placed for south- 
ern observatories than for northern. In 1968 it 
will be favorably placed in both hemispheres and 
probably brighter than the thirteenth magnitude; 
unfavorable will be its short visit of only 12 days 
brighter than 15th magnitude, and rapid motion 
up to 50” per minute. Spencer Jones suggests 
that the latter disadvantage might be overcome 
in part by a variant of the method of interrupted 
exposures.! He recommends also that the color of 
Icarus be determined at an early date. 

(b) For the mass of Mercury, we find Icarus 
closer than Venus only about one-eighth of the 
time, but it can pass Mercury within 8 million 
miles. By chance the next such passage will be 
about 1968 May 1, just before the close approach 
to the earth, It passes within 12 million miles 
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probably every other time; within 20 million 
miles, probably two times out of three. More- 
over, Icarus is sometimes between Mercury and 
the sun, so that the mass of the former is not in 
part confused with the latter as it is in Mercury’s 
perturbations on Venus. The ratio of the periods 
of Icarus and Mercury is 4.645, so that every 
third return of Icarus will find Mercury in some- 
what the same region, with the closest approaches 
at intervals of 14 or 17 returns, the longer inter- 
val being also the 19-year commensurability with 
the earth. Because the coordinates of Mercury 
were not conveniently available we elected to 
postpone the treatment of its perturbations on 
Icarus. 

(c) For the relativity effect on the motion of 
the perihelion of Icarus we have the figure 10705, 
to be compared with Mercury’s well known 
43°03. A recent paper by Gilvarry,? however, 
shows that when the Icarus figure is multiplied 
by certain factors that attempt to evaluate the 
observability of the phenomenon it becomes 
about 14’, as contrasted with Clemence’s esti- 
mate® for Mercury, 3703. These estimates are 
necessarily subjective but suggest that accurate 
observations of Icarus are justifiable. 

2. The ephemerides of Icarus illustrate both 
the special character discussed in Section I and 
the difficulties that its position and speed often 


Figure 1, Positions of Icarus and the sun, 1968 June 6 to 24. 


1953 Aug. Tella or RLOUN OMe, Coat [OU RaN TALS I57 
TABLE I. OBSERVABILITY EPHEMERIDES 

Year Date m p (+34) (—26) a’ cos 6 5’ 

1949 June 29 16.0 0.254 160°E +28° +89° 
uly 7 16.9 0.364 150 E +27 +87 — 34’ =i le 7" 
July 15 17.6 0.480 135 E +26 +86 

1950 July 28 18.9 0.710 80 E +30 +60 
Aug. 7 19.2 0.873 79 E +24 +62 + 69 — 46 
Aug. 17 19.6 1.054 78 E +19 +66 

1951 May 24 20.0 1.501 53 W = +33 
June 13 19.3 Tei 48 W — +29 -+ 100 ar) 
July 3 18.4 1.084 31 W = +12 
Aug. 12 18.6 1.037 37, + 6 +18 
Sept. I 19.4 1283 51 E +12 +31 + gI — 50 
Sept. 21 20.1 1.564 53 E -++10 +34 

1952 June 17 19.9 1.195 78 W = 6 +58 
July 7 19.2 0.875 77 W == % +58 +113 +-- 25 
July 27 18.7 0.667 61 W + 6 +40 

1953 Aug. II 17.8 0.484 114 W 0) +58 
Aug 21 1726 0.411 102 W — +46 +12 et 
Aug. 31 Lay 0.381 83 W — +30 9 
Sept. Io 18.0 0.398 60 W a +20 

1958 May 17 18.0 0.430 70 W +34 +50 
June 6 1720 0.355 119 W +37 +82 Ss ee 
June 26 16.6 0.428 168 W +22 +83 = 
July 16 18.2 0.665 146 E +21 +80 

1959 June II 18.8 0.459 25 E + 2 — 
July 1 17.5 0. 304 74E 45 +50 +141 — 106 
July 21 18.3 0.582 98 E +37 +80 

1968 June 8 17.0 Oni22 32 W +25 — 
June 12 15.2 0.063 50 W +28 — 
June 16 12.9 0.048* canine 1B +88 +32 +240 — 1200 
June 20 13.9 0.100 135 E +52 -+68 
June 24 14.9 0.165 7 Ley +41 +79 


* 1968 minimum distance = 0.042; minimum possible distance = 0.038. 


give to the observer. Table I shows the magni- 
tude, m, the geocentric distance in astronomical 
units, p, the angle of elongation, y, the maximum 
altitude when the sun is 18° or more below the 
horizon for latitudes 34° N. and 26° S., in the 
columns headed (+34) and (—26), respectively, 


TABLE II. 1953 FINDING EPHEMERIS 


1953 a (1950.0) 5 (1950.0) 
Aug. 7.0 12 41™7 +1070 a SAE! ee ae: 
: wh 4I 
aa 2 ae sie Bes 2 —148 
13.0 int hee 508sou eegee 
: i qpitfie® eae —158 
abe 2 ote 122 ays 7 —158 
Pie en So toce | aS SS 254 
aes 4 56.3 eo =GOLA0 Me 
23.0 439.6 143:3 —7EsAGmen aet 
25.0 5 32.5 1 22:9 =73t2T e 1 
27.0 Garg 1094 ely ues 
29.0 Se Li) mia 93 97), 1244 
31.0 8 23.5 +52.4 —70 35 104 
Sept. 2.0 9 6.0 +4as5 —68 23 lee 
4.0 9 39.3 133°3 Stiey od 
6.0 LOM Sed seas — 6290 e202 
8.0 10 2 pee —58 20 qeeelo 
: 5-3 -+16.0 73 -+-212 

10.0 LOeAte 54 48 


and approximate values of the apparent velocities 
a’ cos 6 and 6’ in minutes of arc per day. The 
figures for 1949 through 1952 are shown for com- 
parison with the later ones. Table II shows the 
apparent positions of Icarus for 1953; and Fig- 
ure I, those for 1968. 

3. The observations of Icarus thus far at hand 
and the residuals they show are listed in Table 
III. The orbit was based upon the observations 
of 1949 and 1950, but the representation of 1952 


TABLE III. OBSERVATIONS AND RESIDUALS 
1949 June a (1950.0) 5 (1950.0) Aa(O —C) Ad 


27.23889 16541™08813. —27°51'53"7  —0808 0%O0 
29.23056 16 31 54.65 —28 16 20.5 +0.01 +0.6 


1949 July 

TesisO2me Ol 2401 5.330028 35.05.31 0.08 055 
13.20660 16 03 10.48 —29 26 09.1 +0.03 +0.1 
[4.21701 , 16°02626,90, -—29 28 39.0 ~--0.06, —0.3 
23, TOS83 a 1OL0G OF. 70m 207 47. 25:46) =]-O.03 -—O.1 
24.18750 16 00 06.74 -—29 49 16.5 +0.01 +0.1 
28.20633 16 00. 59.85 —29 56 51.0 +0.02 —0.1 
28.22369. 16 OI 00.01 —29 56 52.4 —0.0I1 -+0.3 
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TABLE III. OBSERVATIONS AND RESIDUALS—Continued 1949 July 28, Mt. Wilson; Baade, Nicholson, Richard- 
son, H, A. C. 1017. 
qosenus: a (950-0) eNNgsoS) Sak It 1950 Aug. 3 to 5, Palomar; Richardson, Nicholson, 
3-19167 13 54 21.03 —I2 OI 25.9 0.00 —I.0 Harrington, H. A. C. 1096. 
3.20972 13 54 25.84 —I2 02 06.5 +0.06 +0.8 1952 June and July, Pretoria; Feast, Lourens, Thack- 
4.19167 13 58 42.56 —I2 39 36.0 +0.06 —I. eray, M. N. Astr. Soc. South Africa 10, 107, 1952. 
BelOrg2 04 02) 56-2813 810105 70m OF Ola Os 
: 1 = 6 42. =— Or : : ; , : 
Be PaO oan) Eno Noes tered oink Sein is so satisfactory as to make differential correc- 
ape done tion unnecessary at least until after 1953 obser- 
BSAA OS4 202 Wan eee ee ae ae vations are available 
23.14514 Ti fa, Life —DOscOnA2 SOjél A : : 
24.13542 I 31 01.7. —I16 15 II —0.6 —3 4. The elements or parameters of the orbit of 
bona Taly Icarus are as follows: 
25.13750 4 29 06.1 — 9 05 03 —0.2 —6 
25.15243 4 29 13.2. — 9 04 35 —0O.I —4 to) =1950 August 7.0 =2433500.5 UT 
ena M=Myotno(t—t)+tAM My=+0.933740 30 
1949 June 27 to July 24, Palomar; Nicholson, Richard- nm =notAn : Mo = 0.015375 5380 
son, Rule, Baade, H. A. C. 1008, 1012, 1017. a, =eP,=a,+ Adz Axo = —0.362753 59 
TABLE IV. PERTURBATIONS IN THE ELEMENTS, UNIT 107° 
Date u - AM Adz Aas Aaz Abz Aby Abz Aa Ae 
1949 
lives — 4 +1093 + 86 +128 + 1 + 49 +23 +167 +17 + 56 
Oct; 11 — 3 + 790 + 48 + 83 —20 + 23 — 3 +120 +75 + 29 
1950 
Jan. 19 2 + 282 =p ly + 45 —17 +3 = 2) + 50 +69 + 16 
Apr. 29 sagt alo es O aS) fo) so. “+4 "25 ao) ag 
Noe, Gf (0) fe) fe) nO) Co) (0) fe) oO oO oO 
Nov. I5 + 1 + 17 = & — 31 —22 — 24 —T4 —=12 7 +14 130 
1951 
Feb. 23 + 2 —- 2 + 24 — 62 —49 — 73 —18 — 46 +43 = 81 
June 3 + 3 —132 7 — 77 —55 —I116 — 9 — 62 +59 —117 
Sept. II + 4 =! 220 + 89 — uf. Ie => RS! fe) —= 50 +30 —i22 
Dec. 20 + 5 = 20%) alos = Oe —48 SST +15 — 48 ape — orn 
1952 
Mar. 29 + 6 102 +117 47 150) Sits neo — 44 ae vl hr 
July 7 Seay = e30 131 35 = 50 — 124 +27 38 +21 — 109 
Oct. 15 + 8 84 +155 — 48 —25 ——T2 1 +25 = ee +18 —II4 
1953 
Jan. 23 +9 — TOF +167 = io —> 6 = iy7/ +43 = AKG — 15 — Or 
May 3 +10 — 102 +156 = Al +6 ROO +52 = 6 SSS) — 68 
Aug. II +11 —- 7 +140 + 9 +18 — 75 +56 — 10 = — 45 
Nov. 19 +12 th 25 +137 bho aia, 3i a AS shee faced + 4 Teer 5 12 
1954 
Feb. 27 +13 =) ie) “133 see +37 =) 72 +59 ad Y, ts — 28 
June 7 +14 — 183 +116 + 2I +48 — 60 +63 — 25 +9 — 10 
Sept. 15 +15 — 209 +103 + 22 +58 — 52 +67 — 137 + 1 + 2 
Dec. 24 +16 — 22] + 99 + 24 +66 — 54 +67 =—EeAG) +30 + 9 
1955 
Apr. 3 eL7 — 400 77 om, +60 = 61 +49 54 72 +. iA 
July 12 +18 — 626 + 69 — 2 +59 — 65 +48 66 +76 fe) 
Oct. ge +19 853 05 — 12 +63 — 68 +49 = BF +65 —- 4 
195 
Jan. 28 +20 — 986 + 55 = By +70 = (oe +48 — 90 +48 een 
May 7 +21 —I140 + 39 — 4I +62 a gt +34 — 102 +72 — 14 
Aug. 15 +22 —1284 + 44 — 57 +49 — 88 +28 —I1I4 +83 — 35 
Nov. 23 +23 — 1503 + 48 — 66 +48 — 97 +29 —125 +80 — 43 
1957 
Mar. 3 +24 — 1683 + 55 = Ge: +50 —I101 32 —129 +68 — 50 
June II +25 — 1763 + 64 — 70 +46 —102 +37 —127 +44 = oe 
Sept. 19 +26 —1791 75 = 983 ao AIG +33 —134 +50 —a7S 
Dec. 2 +27 — 1863 + 88 — 84 +29 — 130 +36 —135 +58 — 83 
195 
Apr. 7 +28 — 1998 + 98 E56 +29 mak, +39 137 +49 = OF 
July 16 +29 —I981 +131 — 62 +37 — 130 +63 — 120 — 5 — 83 
Oct. 24 +30 — 1870 +141 = 45 +41 —124 +76 —113 —23 — 72 
1959 
Feb. I +31 — 1734 +149 me 4 +42 —125 +81 — 109 —15 = 67 
May 12 +32 — 1692 +153 = 320 +44 —128 +84 — 108 — 6 —) 65 
Aug. 20 +33 — 1672 +156 — 25 +51 —123 +89 — 107 —18 — 59 


1953 Aug. 


ay =ePy=a,+ Aa, 
a, =eP,=an+Aa, 


ayy = +0.598282 70 
@20 = +0.440160 94 


oa 
8 
— 

| 


b, =evpO0,=b.+ Ab; = —0.390922 85 
by =evpQy =by0.+ Ab, byo = —0.278029 54 
b. =evVpO0.=b.+Ab, b.o = +0.055733 24 


The differences Aa, etc., with the exception of 
An, are given in Table IV at intervals of 100 days 
and to two fewer places than the foregoing 
constants, together with derived values of Aa 
and Ae defined as follows: 


ao + Aa 
€> + Ae 


Qo = +1.077707 
€>9 = +0.826604 


Il 


é 


From the tabulated quantities we may obtain 
also 


AQ = —59°4 Aa. + 57°9 Aa, — 127°7 Aa, 
+ 60°9 Ab, — 59°4 Ab, + 130°9 Ab, 
At = +13°9 Aaz — 13°5 Aa, + 29°9 Aa, 


+ 39°7 Ab, — 38°7 Aby + 85°4 Ab, 
Aw = —Av — 0.921 AQ 


where 


Av 


+143°3 (d2oAbz + AyoAby + a20Ab-,) 
= —143.3(d2Aadz + dbyoAa, + beoAaz) 


for use with 


Q = Q + AQ Oo = 87.7404 
t= 10 aa At 10 => 22°9788 
w® = wo + Aw wo = 30.9115 


The integration tables from which the quanti- 
ties in Table IV were determined maintain a 
10-day interval throughout most of the calcula- 
tion, with reduction at perihelion passage to 2 
days in 1950 and 1951 and to an even shorter 
interval, because of a close approach to Venus, 
in 1952. Calculations at the reduced interval 
have not yet been undertaken for subsequent 
perihelion passages, but must be before observa- 
tion residuals after 1953 can be regarded as 
entirely consistent with the basic orbit. 

5. Minimum formulae and a sample calcula- 
tion, for problems involving moderate or large 
eccentricities, such as that of Icarus, are as 
follows: 

Given values of M, possibly n, and at least 


five of a2 = ePz, br = eVPQz, ay, by, az, bz, which | 


have been obtained by numerical integration 
from tables, or which are two-body values at the 
start of a step-by-step integration or of successive 
approximations by quadratures. For Icarus, on 
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t = 1952 January 9.0 = 2434020.5 U.T. 
az = —0.36265 b, = —0.39103 
ay = +0.59823 b, = —0.27801 
a, = +0.44011 b, = +0.05569 


M = +2.64564 (radians) 


N= "-1-0.0153755 


Note that 6, may be thought of as having been 
obtained not by integration but from 


A202 + ayby + ab, = 0 (1) 


This equation should be used to obtain preferably 
the 6 that has the largest a as coefficient. I have 
preferred to integrate all of the b’s, however, 
using equation (1) for check. This check is an 
exceptionally good one for finding errors in the 
tables and should be used every ten steps or so 
after enough central differences are available to 
insure full accuracy in the integrals. With x also 
integrated, however, a good check will be found 
in equation (6). 


az +a/,+a/? = e = 0.683091 (2) 


e = 0.82649 
(b.? Fb? + b,7)/E = p =.0.341529 (3) 
Pil ee 10107709) iA) 

a@ = 1.03812 

Genre 11877 


Then, with k = 0.01720209895 for the Gaussian 
constant, the mass of Icarus = 0, and the mass 
of Mercury, m8 = 0.000000 163, included with 
the sun’s mass, 


RVI + mg = km = 0.017202 10036 (5) 
knfat = % = 0.0153759 (6) 


Equation (6) is here used for check, but could be 
a primary source of . Then 


M=E-—esnkE = +2.64564 (7) 


is solved for E by whatever technique the com- 
puter prefers. In this example 


E = +2.86853 
sin EH = +0.26968 


| 


cos H = —0.96295 
ecos H = —0.79587 


a(I —ecos Ff) =r = +1.93539 (8) 
NG 6sin f= 77 = 4 = -+0523138 (9) 
fee eee 50305 (10) 
6 cos A et eae eo 

pee ae a OcAN 122 (11) 
u/r =% = w = +0.11955 qo) 
uw —-wi=e? = 0.683084* (13) 
(ub, — wa,)/e2 = x = +0.71372 (14) 

Fae’ 4 O00 5 

3 = —1.00806 
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(ub, — War) /e? = £ = +0.29887 (15) xy = +4.70316 
y = +0.06266 yy = +1.45432 
2 = —0.11055 zy = +0.50892 
Boba? eg ae y 28 74577 a 16) 1/rs? = 0.008249 28 
xt + yy + 22 = rh = u = +0.23139* (17) m;/Mr;? = 0.000007 8763 
+t 2 = 3 = (18) xy — x = x = +3.989044 
ile 0) (19) yis = +2.94435 
up = xb, + yb, + 2b, * (20) fy = +1.51698 
—W = Xz + yay + 2a, 3 (27) xy? + yi? + 21? = ry? = 26.88606 
G5 =k — Be t (23) ms/Mry> = 0.000006 8488 
bs = we — we : (23) Th the remainder of this tabulation mumenem 


An asterisk indicates a check equation. Where 
no numerical value is given the check was not 
used in the actual calculation. 

The following masses were adopted for the per- 
turbation terms, Mercury being included only as 
augmenting the solar mass: 


Sun + Mercury + Icarus: 


M = 1.000000 163 


Venus: .000002 45 
Earth: .000003 03577 
Mars: .000000 32 
Jupiter: .000954. 79 
Saturn: .000285 58 
Uranus: .000043 73 
Neptune: .OO0005I 78 


The calculation for Jupiter alone based on 
heliocentric coordinates of Jupiter from Plan- 
etary Coordinates, 1940-1960, p. 68, will be indica- 
ted as follows: 


values will be expressed in units of 1078, 


“| a = [#]: = — 972.05 
[y ls = + 871.07 
[@]y = + 638.12 
Lély + [@}e+ [él + [a]; + --- 
= [¢] = — 408 
[gy] =+ 754 
[é] = + 579 
xLé]+ s[9]+ 2Lé] = Lv] = —1998 
2(4L4]+ gly] + 22) = [#4] = — 277 


2ulu | — ld] = [wv] = + 113 


TABLE V. SAMPLE OF INTEGRATION TABLES, UNIT 1078 


D 
Hee [Cn] =[n] 
Dec. 10 — 8498.05 ne 
20 — 8537.24 meee 
30 — 8564.21 pee of 
Jan. 9 — 8581.66 ee 
19 — 8592.22 ‘5 
Date : 3 . ° 
1951-52 2(f] Leal Cf] &Cf] 
Cf] = [M1] 
Dec. 20 ara +3271 ee +128 
Jan. 9 22869 +4269 
Cf] = [a] 
+60756 = 
Dec. 2 461142 iL 386 + 38 cs 39 
] 3 +61566 Spee on Head 3 
ee -+62060 494 
Cf] = [4] 
Dec. 20 cele Saye Te ky 
30 Meera oa’ Fey, Pee ee 
Jan. 9 eet nc aes 2 


Lift = [8] = —1033 
[a ]x — [wl — ulé#] = Laz] = + 494 
ri = =f 363 
a — ee 
[w |x — [w]e — we] = [b. |= — fren 
[o,] = a 
[6.]=+ 131 
[rn] 6[n] &[n] [2] 
+12.22 cea a +.24 ee 
—2.70 —.I 
a Aen 2203 +.07 é 
al =f] Cf) een ef] asf] 
Cf] a [a,] 
—16 — 27903 a — 00 + 6 
= ib 78 = ap 7 
e 2 od 
— 27906 f 
Lf] = (6.] 
— 8 — 64685 117 + + 8 
“7 3 on 
Cf] = re 
+10 — 27856 ny — 72 + 4 
es hes, +25 
are Eetiee +151 7 47) ta ae 
—27376 9 113! 
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(az bz] + a,[by] + a.[b.])/e differences determined, so that the ends of the 
= [bo ]vVp = + 438 tables appear as in Table V. 
— (b,[ dz] + b,[a,] + 0.[ 4. ])/e The table for [b, ] is not included among the 


ca [ove = 4+ 437% foregoing because it is not strictly necessary in a 


’ ; ui yc minimum calculation, but it has always been 
(Lé]Vp — 2[w))/Va = [M]= Wade carried as a check, and I| strongly recommend 


~analé|=[%]=+ 6.89 that it should be so carried 
~ (x[a2] + yLay] + 2L4.]) = [vb] = i 112" Correction to M and the starting elements are 
(x[b2] + yLby] + 2[6.] + ulw])/p obtained from these tables by means of the fol- 
a ‘a = [u : lowing formulae: 
2(éLde]+ ilay]+ Ala] + ela) /u ; : 
= [a] ; AM. = kn { { [nee bn { - [at at 2 
~{w[u] 4 Laul a |(p 7 r)} Jase z m F [7 | + Rm 2 [ ] ( 4) 
= [M]= +426 inary) 
Soe ores Afi; = Rm lef Eben 09h e Cena, CfCe 25) 
The variations [7], [M1], and at least five of : 


[a.], [ay], [a], [b.], [6,], [6] are now entered where, for step-by-step integration, when central 
into the tables, and the sums and backwards _ differences are not available, 


Rin a [7 |d#2 = wk 


Rm ‘A ate = WRm {2c + + 3[f]i- 1 ts ray > fla + = ah Cf ]i-2 Hee 2 sft 4 of, 
f= Mn, C705, CLC 27) 


Me eae sas 5 Die ere flea + Mites + 5 Seat | (26) 


In the present example only one, two, or three 6. Comparisons between the method of varia- 
terms of these expressions are significant. With tion of parameters and the methods of Cowell 
the initial elements for tp = 2433500.5 UT, as and Encke, and between the formulae of Section 
shown in Section 4, and 5 and older formulae for the variation of param- 

iis AL gts AEH AME (28) eters will illuminate special perturbation theory 


and its practical application to objects such as 
f= fot Af, f=, az, be, etc. (29) Icarus. 


the example yields, on 1952 January 19.0 = Cowell’s method, as is well known, involves 
2434030.5 UT, the integration of derivatives that are larger and 
AM = —0.000183 M = 2.79941 that have larger differences than those of Encke’s 
Pe 6 00000001 nm =  0.0153755 method or the method of variation of parameters. 
Aa, = +0.000107 a, = —0.36265 The contrast becomes more pronounced, if the 
Aa, = —0.000053 dy = +0.59823 orbit is eccentric, as the object approaches peri- 
Aa, = —0.000048 Gg =o A40l8 helion. In these circumstances Cowell’s method 
Ab, = —0.0001II b, = +0.39103 is customarily replaced by one of the other two, 
b, = —0.27801 but it is profitable, nevertheless, to consider some 

Ab, = —0.000047 b, = +0.05569 of the differences that would appear if we at- 


tempted to bring an integration of Icarus to 
and the calculations are ready for the next step. perihelion without diminishing a 10-day interval. 
It will be noted that the foregoing derived values, Table VI supplies the second derivative of x and 
except WM, are identical with those used at the its differences as they would appear in such an 
start of the step illustrated. It was necessary in attempt by Cowell’s method, together with the 
fact to use the integration formulae only for derivative of one of the parameters and the differ- 
every fourth or fifth step of the calculation. The ences of this derivative as they appear in the 
interval of the tables was determined by the run method of variation of parameters. 

of the differences and not by interpolation or It will be noted at a glance that the sizes of 
integration. the derivatives and differences favor the latter 
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method. Draper found it impossible to carry 
Icarus by perihelion with an interval one-fourth 
as great, in the work upon which the first calcu- 
lations of the present paper were based.’ 

My purpose in reproducing Table VI, how- 
ever, is to show not so much the contrast between 
the sizes of the integrands and their differences, 


TABLE VI. COWELL VS 
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as the similarity of the oscillation patterns of the 
two sets of differences. The rapid change in the 
position of Icarus at perihelion, 180° in about 
10 days, is clearly reflected in the perturbations 
as well as in the principal or solar terms of the 
accelerations. The differences of the variations of 
the parameters, though small, are yet erratic 


- VARIATION OF PARAMETERS 


Rm? dt? j & 53 54 
June £3 0)— 1+0053907.55 ae —On2 704 Ie en — 0.3874407 
23815575 LOO 2m ee —EO.c25354 1 seme Key — i007 7 b2 71 
July; 132) 482 -97025597 een 563, 9°0392397 444 7493799 36.9632855 
13, — 7.4046421 unite ses 3 o> Le LOLA Oma rae A — 109.8118522 
23 +19.8711418 7-2757°39 — 43 .3523321 75: 4323 +131 .7093642 
I day 
ie Ge 6 62 53 64 
June13 =—0.00000230 — 69 + 35 
23 —0.0000022I a 7 — 55 es ap Wshe) 
July 3  —0.00000267 - . + 95 ae + 6 
13 —0.00000218 o a +251 wi oa a7 
23, + +0.00000082 3 —329 5 + DI51 


enough in their localized oscillation at the selected 
10-day interval to cause uncertainty in subse- 
quent values of the integrated parameters, and 
we are forced to a 2-day interval to get by peri- 
helion. It is evident, however, that the uncer- 
tainty is too small to affect the calculated values 
of the variations for some distance beyond peri- 
helion. We are thus able to carry well past with 
the 10-day interval, compute all at once the 
values of the parameters at the shorter interval, 
then integrate by quadratures to obtain the nec- 
essary correction to the second sum or integral 
in the 10-day tables. - 

The amount of labor necessary to find correc- 
tions of about 2 units in the seventh decimal 
place, due to the rapid motion of the perturbed 
body rather than to its proximity to disturbing 
bodies, nevertheless indicates that it is urgent to 
find other solutions to the difficulty. Three 
avenues deserving of exploration are: (1) In the 
mathematics of finite differences to investigate 
means to correct for the effect of the oscillating 
differences of a near-singularity or the extent to 
which they can be tolerated. (2) On the astro- 
nomical side to attempt to extend the variation 
of parameters to include a part of the perturba- 
tions as well as the solar terms of the equations 
of motion. (3) The substitution of another inde- 
pendent variable for t, to which possibility I shall 
devote some discussion in the next section. 

A 20-day or 40-day interval would have been 
possible in the neighborhood of aphelion, but the 
10-day interval was continued because of the 
shortness of the period of Icarus (409 days). 
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Two further comparisons with Cowell’smethod: 
(1) The computers found the labor of computing 
the variations of the parameters from the per- 
turbations small in comparison with that of 
computing the perturbations themselves from 
the coordinates of the seven disturbing planets. 
(2) First-order perturbations cannot be obtained 
by one double-summation “‘quadrature”’ in Cow- 
ell’s method ; at least the principal or solar terms 
must be determined step-by-step or by succes- 
sive approximations. In the variations of the 
parameters, on the other hand, these terms are 
eliminated, so that first-order ‘‘quadrature”’ inte- 
grations without successive approximations are 
possible. | 

The foregoing two comments apply also to a 
comparison with Encke’s method, whose periodic 
rectifications and gradually increasing (fgx — &) 
terms should be balanced against the calculation 
of the variations of the parameters. Otherwise, 
in step-by-step integration, the problems of dif- 
ferencing and intervals will be much alike in the 
two methods. In first-order integrations, how- 
ever, Encke’s method is in the same relatively 
disadvantageous position as Cowell’s method. 

The first and most obvious comparison between 
the formulae of Section 5 and the older formulae, 
as shown by Watson,® Oppolzer,® Bauschinger,’ 
and others, is that the former refer to trigono- 
metrical tables only for functions of the eccentric 
anomaly while the latter refer also for v, w or u, 1, 
and, for each disturbing planet, 6 and / — Q. 
The integration of az, b,, etc., moreover, elimi- 
nates much of the calculation involved in the 
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combinations of sines and cosines into the com- 
ponents of orientation vectors. Perhaps even 
more important are the checks supplied by the 
integration of x and the third 6; no such checks 
have been used or are readily available with the 
older formulae. 

In problems involving a resisting medium, of 


iH Ba AS DROWN OMECAL. JO URINATE 


which there may be more in the future, the new 


formulae are by chance at an advantage in that 
they involve the calculation of the components 
of the velocity, @, y, 2, rather than the compo- 
nents of the transverse vector (in the orbit plane 
at right angles to the radius vector). 

The way in which each set of formulae meets 
the special requirements of very small eccentrici- 
ties is a matter of some interest, though better 
reserved to another paper than discussed at 
length here. | shall remark only that both formu- 
lations replace the mean anomaly, M, by the 
mean longitude, L, because [M] has ane in the 
denominator, effectively first power only, and 
LL ] does not. In the formulae of Section 5, more- 
over, we must replace the equations for [6, ], [b, ], 


[6.] by 
[é.] = yL2] — 2Ly] 


rad aa) — [2] (30) 
Lé.] = x[y] — 94] 
where c, = VpR;, etc., are free from the e in 


= eVpQ;, etc., so that the components of the 
orientation vectors may be more accurately de- 
termined. The uncertainty remaining ina, = eP,, 
etc., has its counterpart in the older formulation 
in the determination of 2; the position of peri- 
helion is thus made uncertain, but not the posi- 
tion of the object under investigation. In this 
formulation I have been considerably anticipated 
by Milankovitch.$ 

7. Modifications or additions in or to the mini- 
mum formulae of Section 5 may be sought to 
avoid certain disadvantages, but will in turn 
introduce certain others. Before discussing pos- 
sible changes I list both the disadvantages present 
in the minimum formulae and those introduced 
by changes: 


a) No check upon the integration of M. It 
will be noted that equations (1) and (6) pro- 
vide checks for @z, dy, @z, bz, by, bz, and n, but that 
there is no final check, other than recalculation, 
on the integration of [M] or even of [7] by (27) 
and (26), respectively. 

b) Solution of Kepler’s equation. I should note 
that the computers to whom I give credit at the 
end of this paper found Kepler’s equation un- 
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objectionable, especially as compared with the 
devices proposed to avoid it. 

c) Changes of interval required by perihelion 
passage rather than perturbations (cf. Section 6). 

d) Required interval even shorter than that 
of the minimum formulae. 

e) Interpolation of planetary coordinates when 
the arguments for the tables do not coincide with 
the dates in Planetary Coordinates. 

f) Said interpolation not possible in advance. 

g) Reinterpolation necessary for each approxi- 
mation. 

h) Accurate evaluation of ” required before 
the tables of / or t can be extended. 

Table VII indicates how these disadvantages 
are found in the several devices that have come 
to my attention. Brief notes are given describing 
these devices, but a more complete analysis of 
the evaluation is left to the reader. 


TABLE VII. COMPARISON OF ADVANTAGES AND DISAD- 
VANTAGES IN VARIOUS INTEGRATION PROCEDURES 


d e-f ‘¢ 


~ 
i i oe 


. Minimum formulae, as above. 

. Integrate Mp» or T. 
I+ II. 

. Integrate £ or other variable. 
I + III. 

. M as independent variable. 
II + IV. 

. Eas independent variable. 

. € as independent variable. 


i i i i i <i <i e) 


mw MMM 

mM mM 

mM mM 
tal 


NOTES ON TABLE VII 


IJ. The integration of Mp or 7, treated as parameters, 
is based upon the equations, 


[Mo] = [M] - [n]¢ — b), (31) 
(T] = {l#]¢ — T) — [M]}/n. (32) 


The latter of these equations I included in my first dis- 
cussion of the method of variation of parameters.® The use 
of either is compared to the usual procedure, as set forth 
in Section 6 above, in my more recent paper.!2 We made 
use of one or the other of these procedures to provide a 
check upon M throughout most of our work on Icarus. 

III. The integration of the eccentric or of the true 
anomaly for Icarus, to avoid the solution of Kepler’s equa- 
tion, cannot be carried on with an interval larger than 
can be used with Cowell’s method. Experiments with 
N = E—esin N, and with other such variables whose 
differences might be expected to be smaller than those of 
E, were disappointing. 

V. The use of M as independent variable has been dis- 
cussed recently by Merton." The perturbative variations 
are to be multiplied by, and ¢ is to be integrated from 


Jaen, OH) 
{1 = knLM ] am 


as Merton formulated it, or, more simply and avoiding the 
possible need for a second approximation, by 


A= (nt bal] 


This gee: avoids the solution of Kepler’s equation by 
tabulating certain functions for integral values of M and e, 


(33) 


n 


(34) 
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as proposed originally by Crommelin for first-order quad- 
ratures.”” 

V. The use of E as independent variable requires that 
the perturbative variations be multiplied by and ¢ be 
integrated from 


r 
a 


{n + kml M]+kn[e] sin E}-. 


The differences are well behaved and the interval of ¢ is 
automatically shortened around perihelion. 

VI. If we designate by « the eccentric anomaly in a fixed 
reference orbit, let us say the initial osculating one, ¢ is 
obtained without integration or approximation from 


(35) 


e— éosine — My 


t—ft= ae (36) 
and the perturbative accelerations are multiplied by 
dt 1 — eocose 
dca pier iawn: (37) 


Both M and E would be integrated to check through 
Kepler’s equation. For comparison (orbit of Icarus): If the 
interval of « is 10°, there are 36 steps per revolution and 
the interval of ¢ varies from 1.96 days at perihelion to 20.7 
days at aphelion. In our calculations by the ‘‘minimum 
formulae” there were 41 ten-day steps per revolution plus 
about 25 two-day steps to carry the integrals by perihelion, 
a total of 66 steps. Since the 25 two-day steps also required 
interpolation, though of a somewhat simpler character 
than that considered under disadvantage e), it seems that 
the use of e as independent variable should be seriously 
considered at least for objects whosc orbits are as eccentric 
as that of Icarus. In first-order quadratures e is synony- 
mous with £ and the proposed device, requiring no tables, 
appears to be considerably simpler than Crommelin’s.” 
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